Rigorous analysis of bistable memory in silica toroid microcavity 
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1. Introduction 

Optical bistability is a fundamental physical phenomenon 
where it is possible for certain devices to have two stable 
transmission states. It occurs when the refractive index 
or the absorption of the nonlinear medium in an opti- 
cal cavity is dependent on the light intensity. Optical 
bistable devices are considered to be an important build- 
ing blocks in all-optical signal processing for such compo- 
nents as optical memories and optical flip-flops [l(, and 
this phenomenon was extensively studied in 1970s to 80s 
0. However, at that time, the size of the cavity and 
the operating energy were too large to be considered for 
practical applications. 

Recent progresses in achieving higher quality factors 
(Q) in ultrasmall microcavities on-chip has refo- 

cused attention on optical bistability [7| , due to the pos- 
sibility of achieving denser integration and lower energy 
consumption. Since the photon density in a cavity scales 
with Q/V, where V is the mode volume, a high Q cavity 
with a small V enables us to use various nonlinearities 
at extremely low input powers; hence it allows optical 
bistability at an ultralow power. 

Soljacic et al. demonstrated numerically that optical 
bistability based on Kerr nonlinearity is possible by us- 
ing a two-dimensional photonic crystal nanocavity at an 
ultralow driving power of 133 mW [7]. Various experi- 
ments have already been reported in silica microspheres 
, silicon photonic crystals d, [To[ , and silicon micror- 
ing resonators [III, but by using the thermo-optic (TO) 
effect. Since the TO effect is accompanied by thermal 
accumulation, the response is relatively slow. To achieve 
faster speed, the carrier-plasma effect has been utilized. 
This effect is the result of carrier generation [ll| fl3j . 
The latest research has reported on 4-bit optical random- 
access memory operation based on carrier nonlinearity, at 
a power consumption of only 30 nW, by using InGaAsP 
photonic crystal nanocavities [l4j . The key to achieving 
the low power consumption is the smallness of the cavity 
V. 



As introduced above, the required power for the oper- 
ation of optical bistability has been significantly reduced 
in recent years due to the high Q/V. However, all of 
the demonstrations use either TO or carrier-plasma ef- 
fects to drive the bistability whose nonlinearities are ac- 
companied by photon absorption. Recent advances on 
linear and nonlinear studies of microcavities and chip- 
based waveguide devices have opened possibility of their 
use not only for classical all-optical processing but also 
for loss sensitiv e appli cations such as quantum informa- 
tion processing |15l4l8l |. With those applications in mind, 
we need to reduce significantly the power loss (consump- 
tion) of the bistable system. The use of the optical Kerr 
effect is the ultimate goal because it does not absorb pho- 
tons. In addition, the use of an ultra-high Q allows us to 
reduce the power scattering loss of the signal light from 
the system, and this will support all-optical information 
processing for loss-sensitive applications. 

Although there have been several attempts to employ 
the optical Kerr effect by using large-bandgap materials 
such as AlGaAs and chalcogenide glasses [20| , it still 
seems to be difficult to obtain Kerr bistability without 
suffering from the carrier effect or TO effect jg, [2l|, [22j ■ 
On the other hand, silica has been an excellent mate- 
rial with which to study of various aspects of based 
physics [23|], because of the large bandgap that can sup- 
press carrier generation. Therefore silica microcavities 
have the potential to achieve a Kerr bistable memory 
that consumes very little energy. 

Among various silica microcavities, the silica toroid 
microcavity [5J| has an ultra-high Q and is capable of 
integration on a chip. As discussed above, a high Q cavity 
is attractive for both achieving optical Kerr bistability 
and for low-loss applications. 

In this paper, we demonstrate numerically that an op- 
tical bistable memory based on the optical Kerr effect 
is possible by controlling the coupling between the cav- 
ity and the tapered fiber. Our model consists of cou- 
pled mode theory (CMT) and the finite element method 
(FEM). 
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The paper is organized as follows. First in $2j we pro- 
vide a simple picture of the physics that we are going to 
employ. In $3] we describe a numerical simulation model 
that combines CMT and FEM. In §4 we show the calcula- 
tion result and in §5 we discuss the energy consumption. 
Finally, we finish with a conclusion. 

2. Simple model 

As described in the previous section, it is still difficult to 
use the optical Kerr effect in ultrahigh-Q silica microcavi- 
ties, even though the material has a large bandgap. This 
is because of the small light absorption at the surface, 
which is caused by a thin water layer 0, [32j. There- 
fore, we need to analyze this carefully by using rigorous 
modeling and applying realistic physical parameters to 
reveal the conditions required for obtaining Kerr bista- 
bility. However, before undertaking a rigorous analysis 
we start with a simple model to gain an intuitive under- 
standing of the strategy for obtaining Kerr bistability. In 
this section we pursue an analytical discussion about how 
we can obtain the Kerr effect without exhibiting consid- 
erable TO effect. 

First we derive the relationship between the wave- 
length shift and the input energy required for Kerr non- 
linearity. The energy U of an electromagnetic wave in a 
dielectric medium is given by 24 \, 



U 



1 



eEtV, 



(1) 



where e, E$ and V are the dielectric constant, the electric 
field amplitude and the mode volume, respectively. Using 
Eq. fTJ and / = ve\E \ 2 /2, where v is the light speed of 
the medium, we obtain the relation between the power 
density / and the optical energy U as, 



2c U 



(2) 



where c is the velocity of light and no is the refractive 
index of the cavity medium. Then, the refractive index 
change caused by the Kerr effect is given by, 



An K err = n 2 I = 



2n 2 cU 



(3) 



where n 2 is the nonlinear refractive index. Equation ([3j) 
allows us to calculate the refractive index change as a 
function of the energy in the cavity. 

Next we discuss the TO nonlinearity. The refractive 
index change caused by the TO effect Anxo is given by, 



An T o = n £T. 



(4) 



Where £ = (l/n)(dn/dT) is the TO coefficient. The 
energy required to increase the temperature of 1 K for 
volume V is given by CpV, where C is the heat capacity 
and p is the density of the material. Then, using Eq. (0} , 
we obtain 



An T o 



np£, 1 
Cp V 



U&hs, 



(5) 



where £/ a b s is the energy absorbed by the material. Equa- 
tions ([3]) and ([5]) describe the refractive index change at 
a given energy for the Kerr and TO effects respectively. 
For the Kerr effect to be larger than the TO effect the 
following condition is required. 



An K crr 

An T o 



2n 2 cpC U 
n?£ C/ abs 



> 1 



(6) 



To gain a simple understanding, we assume a steady- 
state model. The light energy in the cavity is constant 
and generates heat at a constant rate. Thus, U(i) and 
U & bs(t) are expressed as, 



U(t) = U 
d£W*) = L 



u. 



abs 



1 



U(t), 



(7) 



where Uq, T^is and r a b s are the steady energy in the cav- 
ity, the thermal relaxation time and the thermal gen- 
eration rate caused by photon absorption, respectively. 
When we neglect the thermal diffusion (r^ = 0), Eq. (|6]) 
is simply expressed as, 



2n 2 cpC 

* ^ o~Z Tabs 

n ot 



(8) 



This provides direct view as to how we should design the 
cavity system in order to obtain the Kerr effect without 
the TO effect being too great. When we use the following 
parameters for Si0 2 , £ = 5.2 x 1CT 6 K _x , n 2 = 3.67 x 
10 -20 m 2/ w n = 1>47) C = 741 x 10 -l J/( g -K), and 

p = 2.65 g/cm 3 , Eq. (JSj> gives the following condition, 



t < 3.84r : 



abs 



(9) 



This equation provides a simple understanding of how 
we can achieve Kerr nonlinearity in Si02 microcavities. 
Although the Kerr effect governs the refractive index 
change (AnKerr > Anxo) at an early stage of the op- 
eration, the TO effect becomes dominant after a period 
of time given by Eq. ©. When we employ a r a b s of 329 ns 
(the selection of this value is discussed in detail in §3 C|) . 
Eq. © gives t < 1.26 ps, which is the duration time, 
for which AnKerr is larger than Anxo- Note that even 
if we consider the thermal relaxation time in our model 
as, Tdiff = 8 ps (obtained from Fig. [3]), this duration time 
does not change greatly and is about 1.35 ps. 

To achieve Kerr nonlinearity we must complete our op- 
eration before this time, namely the time that the light is 
absorbed by the material. This simple picture is straight- 
forward to understand. If we are to obtain a faster opera- 
tion speed we require a smaller rtot because it determines 
the rise and fall time of the light energy in the cavity. r to t 
is given by r^ 1 = ^ + r"^ + t~^ s , where r loss and r coup 
are the photon lifetimes denned by the loss rate that do 
not contribute to the generation of heat and coupling to 
the waveguides. Since r a b s and ti oss are mainly deter- 
mined by the material and the structure that we use, the 
only parameter that we can control is r coup . This dis- 
cussion suggests that optical Kerr operation is possible 
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by controlling the coupling between the cavity and the 
waveguides, because it enables us to release light into the 
waveguides before it is absorbed by the material. 

In the following section, we perform a rigorous anal- 
ysis showing that Kerr operation is indeed possible by 
changing T coup . 

3. Rigorous modeling of the optical Kerr effect 
and thermo-optic effect in a toroid microcavity 

A. CMT in whispering gallery mode resonator 

First we describe our master equation based on the cou- 
pled mode theory in a whispering gallery mode (WGM) 
resonator to obtain the linear and nonlinear transmit- 
tance. The structure is shown in Fig. Ufa). Because a 
two-port system (a side coupled cavity with one waveg- 
uide) makes the bistable operation difficult to observe, we 
focus on a side-coupled four-port system [25] throughout 
this paper. It consists of a toroid (ring) cavity and two 
waveguides for input and output light, a detailed discus- 
sion on the comparison between side-coupled two- and 
four- port systems will be provided elsewhere [2(1] . Briefly, 
as discussed in $2] we need to make the coupling large in 
order to prevent heat accumulating in the system. How- 
ever, the transmittance spectrum of a two-port system is 
shallower in an over-coupled configuration, which makes 
the switching contrast of the output light very low and 
difficult to distinguish. Hence, optical-bistable switching 
with high contrast is difficult to achieve in a two-port 
system in the presence of the TO effect. 
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Fig. 1. (Color online) (a) Cavity structure used for nu- 
merical analysis, (b) Cross-section of the mode intensity 
profile of a toroid microcavity. [27| 



The mode amplitude a in the cavity is given as [28[ , 

1 1 



da 
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2 \ Tabs ^loss ^coupl 7~( 



coup2 



l 



7"coupl 



exp(j9)sh 



(10) 



where uq is the resonant frequency of the cavity. The in- 
put wave s; n excites the counter clock-wise (CCW) mode 
in the cavity. We assume an ideal cavity where there is 
no coupling between the CW and CCW modes. r coup i 
and T coup2 are photon lifetimes determined by the cou- 
pling with the lower and upper waveguides, respectively. 
9 is the relative phase between the mode amplitude in 
the cavity and the optical wave in the lower waveguide 
and is given as, 



J7r-//„(i? + /i [^-J 



(11) 



Here R, r, A and Ao are the major and minor radiuses 
of the cavity (shown in Fig. [lib)), the input wavelength 
and the resonant wavelength of the cavity, respectively. 
Equation (|lip shows that the phase between a and Sj n be- 
comes unmatched on off-resonance. Output waves s ou ti 
and s ut2 are given as, 




exp(-j6q)a 



exp {-jf3 2 d) 



7coup2 



(12) 



(13) 



where fix and /?2 are the propagation constants of the 
lower and upper waveguides and d is the waveguide 
length. Note that the relative phase between the mode 
amplitude in the cavity a and the upper waveguide s- m 
is always zero because there is no incident wave in the 
upper waveguide. 

When we use a slowly varying envelope approximation, 
i.e., 



a(t) = A(t) exp (jut) 
s-m(t) = Si n (t)eyq>(jut), 

we can rewrite Eq. (fTU)) as, 



(14) 
(15) 



dA(t) 
dt 
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2 V T abs 
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1 

'Hoss 



1 
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+ 

^coup 1 T" CO up2 



1 



7"coupl 



exp(j6)S- m (t). 



A(t) 
(16) 



where A(t), S- m (t) and u = 27rc/(n A) are the envelopes 
of the cavity mode and the waveguide mode and the fre- 
quency of the input wave, respectively. Equation (|16p is 
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the master equation of the linear system. By using this 
equation, we now can calculate the energy in the cavity 
and the output power at an arbitrary time. 

B. Modeling the nonlinearities 

To describe the nonlinear effects in our model, we take ac- 
count of the nonlinear refractive index modulation caused 
by the Kerr effect (ArtKcrr) and TO effect (An-ro) hi the 
master equation (Note that the carrier-plasma effect is 
negligible in silica due to its large bandgap). The nonlin- 
earities in an optical cavity result in a shift in the reso- 
nant wavelength because the optical path length changes. 
Thus, the shift of the resonant wavelength of a cavity is 
given as, 



5X(t) 



An(t) 



An 



no 



(17) 



where An(t) is the effective nonlinear refractive index 
change of the cavity. By substituting Eq. (jTT)) into 
Eqs. (IT^|) and (TTTj) . we obtain 



dA(t) 
dt 



2ttc 



1 



no + Ara(t) \Xq 
- (— + — 

2 \ r abs 71 



5X(t) 
1 
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"7~cou.pl 



■exp(j0)S in (t). 



A(t) 



(18) 



9 = 4n 2 (n + An(t)){R + r) 



A + SX(t) 



X ) • ^ 



These are the master equations that we used in our 
model, which take the nonlinearities into account. 

Next, we describe how we calculated the nonlinear re- 
fractive index change Anxcrr and An-ro ■ We can directly 
calculate AnKcn- from Eq. ([3]). Taking the spatial depen- 
dency into account, we obtain, 



An KclI (x,y,t) = n 2 I(x,y,t) 

2n 2 c - 

= Up{x,y,t), 

n 



(20) 



where U p (x, y, t) is the energy density distribution of the 
cavity mode in x, y cross-sectional coordinates. It is given 
as, 



U p (x,y,t) 



2ttR 



(x,y), 



(21) 



where U p — \A(t)\ 2 is the energy of the light stored in 
the cavity and I(x,y) is the normalized cross-sectional 
power density distribution of the whispering gallery mode 
obtained by FEM (see Ref. [13] )• It is normalized as 
JJ I(x,y)dxdy = 1, and the profile is shown in Fig. [TJb). 

The refractive index change An-ro j which is induced 
by the TO effect, is described as, 



An TO (x,y,t) = n£[T(x,y,t) - 300 [K]; 



(22) 



The cross-sectional temperature distribution is calcu- 
lated by using 2D-FEM (COMSOL Multiphysics). By 
setting the heat source in the dielectric cavity as, 



Q'(x,y,t) = 



'abs^P 



U p (x,y,t), in cavity, 



in air, 



(23) 



where r~ bs is the thermal generation rate caused by 
the material absorption, we can obtain the temperature 
T(x, y, t) at any time and at any position by performing 
an FEM calculation. 

Finally, we obtain the effective nonlinear refractive in- 
dex change An(t) as, 



An(t) 



JJ [An T o(z, y, t) + An Kc „(x, y, t)} I(x, y)dxdy 
JJ I(x, y)dxdy 



(24) 



Now by solving Eq. flT8"|) sequentially using Eqs. (fTTj) 
(l24l) . we can obtain the light energy in the cavity U p (t) 
and the output powers P D uti = |Souti| 2 and P ou t2 = 



C. Determining the absorption and the photon lifetimes 
The total photon lifetime T to t is defined as, 



r tot 



cont 
x -1 



rad 



coupl 



coup2/ 



(25) 



where 7" ma t, T wa ^ er) T CO nt7 ^rad; ^suxf: ^coupi and 7" COU p2 are 
photon lifetimes determined by the absorption of the ma- 
terial, water absorption that is usually present on the 
surface of the cavity, absorption caused by surface con- 
tamination, radiation loss, scattering loss, coupling to 
the lower waveguide and coupling to the upper waveg- 
uide, respectively. Since this expression is very compli- 
cated, we define the photon lifetime that is related to the 
absorption as, = r m l t + T~* tel + r~* t , losses to the 
outside of the cavity Tj~ ss = r r ~ ^ + r s ~ J f , and the couplings 
Tcoupi,2- By using these photon lifetimes, we can simplify 
our analysis, because r a b s contributes on the generation 
of heat, but the other factors (ri oss and r^upi^) do not. 

Here we describe the photon lifetimes that we used in 
our analysis. The material absorption of silica at tele- 
com wavelength is usually very small (a — 0.2 dB/km 
[2^]). but, it is known that silica toroid microcavities have 
much larger absorption due to the water layer and con- 
tamination on their surfaces. And the Q factor is limited 
by these absorptions in practice Q. Thus, we decided to 
consider two cases, which we call an ideal case and a re- 
alistic (worst) case. In an ideal case, we assume that the 
absorption loss is determined by only the material (this 
means that the absorption loss is extremely small) and 
the experimental Q is limited by the losses to the outside 
of the cavity t^ ss . On the other hand, we assume that 
the experimental Q is limited by the absorption loss T~ bs 
in a realistic case. This implies the worst case in terms of 
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thermal accumulation, because a large part of the inci- 
dent light eventually turns into heat. Any result should 
be better than that obtained with this latter condition. 

To consider those two cases, first we fix the intrinsic 
photon lifetime r mt = (t ~C + r io SS ) _1 as ^29 ns (corre- 
sponding to Q; nt = 4x 10° |30l| . where Qi nt is the intrinsic 
Q), since this is the record largest experimental T- mt . In 
the ideal case we disregard losses than the intrinsic ma- 
terial absorption; so we use r a ,h s = 164 /is (corresponding 
to Q 



abs 



; mat 
loss 



» r. 



2 x 10 

-l 



li 



3 li p and t\ 



abs ' 



; 330 ns; 

On the other hand, we use 



Le - T int 

Tabs = 329 ns (corresponding to Q a b s = 4 x 10 8 [30|, |32j) 
and Tioss = co for the realistic case. This Q a bs value is 
the same as the highest experimental Q[ n t, for the reason 



» r. 



los; 



discussed above; i.e. r;~ f 

Although Ti n t is determined by the material and struc- 
ture, we can control r coup i and r coup 2 by adjusting the 
distance between the fiber and the cavity [33| ■ With this 
in mind, we conducted numerical simulations for various 
Tcoup2 values. Note that r coup i is controlled to satisfy 
2 )~ 1 and thus achieve critical cou- 



' coupl 



coup2/ 



pling [28, 34j,[3jf between the cavity and the lower waveg- 
uide. In this condition, the power transmittance through 
the lower waveguide T out i = P uti/Pin decreases to zero 
on resonance and high contrast can be obtained between 
two output states [26|. 

4. Numerical calculations 

A. An ideal case: Small material absorption 

In this section, we show that Kerr bistable memory is 
easily feasible, without careful adjustment of r coup , if only 
the inherent material absorption of silica is present. 

First, we input a rectangular pulse to investigate the 
refractive index changes AnKerr and Atito ■ The result is 
shown in Fig.[2][a), where r coup 2 is equal to Tint- It shows 
that AriKcrr is always larger than AriTOj which tells us 
that the influence of the absorption induced thermal gen- 
eration is nearly negligible. Hence the Kerr effect is easily 
obtained without it suffering from the TO effect. 

Next, we input a triangular pulse to investigate the re- 
lationship between the input and output of the system. 
To allow us to charge and discharge the cavity gradually, 
we set the pulse rising/falling rate of the triangular in- 
puts at dPi n /dt = 62.5 nW//is. Figure [2jb) and (c) are 
plotted from the input-output response of a triangular 
input of different detuning values S. The lower coupling 
photon lifetime r coup i is set equal to (r^ 1 + T c ~J lp2 ) _1 to 
achieve critical coupling. When 5 is greater than 20 fm, 
clear hysteresis is observed, which is direct evidence of 
optical bistability. Optical bistability is observed in the 
longer side of the wavelength detuning, because the Kerr 
effect increases the refractive index (shifts the resonance 
toward the longer wavelength). Note that we obtained a 
large contrast between the two bistable states in Fig. &b) 
because the cavity transmittance P ou ti falls to zero on 
resonance in the critical coupling [35J. 

Finally, we performed optical memory operations as 



shown in Fig.[2](d). Again, we set T COU p2 equal to rj n t and 
the detuning S equal to 27 fm. The solid line is an input 
with a drive power P^" vc of 3.2 /iW. The peak power 
of the 0.8-/xs square set pulse is P^ 1 — 5 fiW. To reset 
the system, we reduce the input power to P*° sct = 2 fiW 
for a duration of 0.8 /is, which we call a reset pulse. The 
duration of the negative reset pulse must be longer than 
the discharging time of the cavity, which is equal to T to t ] 
otherwise the cavity does not reset. Set and reset pulses 
are inputted at t = 8 and 20.8 [is. P ou t2 (shown as the 
broken line) rises to high (ON) state when the set pulse 
is inputted. It keeps the ON state until the reset pulse 
is injected. After the reset pulse has been entered, P ou t2 
drops to low (OFF) state and holds this state. P ou ti 
(indicated by the dotted line in Fig. HJd)) shows the in- 
verse behavior of P ou t2- Figure [^d) clearly shows optical 
memory operation, which is based on Kerr nonlinearity. 
Although the "memory holding time" demonstrated with 
this calculation (shown in Fig.[2][d)) is about 30 /xs, it can 
be much larger, since P ou ti and P ou t2 exhibit almost an 
plateau response due to the small material absorption. 

B. A realistic case: Large material absorption 



As shown in §4A] the realization of a Kerr bistable mem- 
ory is feasible without it suffering from the TO effect, 
even when the coupling is not large, if only the inherent 
absorption of silica is present in the cavity. In reality, 
however, other sources of absorption occur in a microcav- 
ity, such as surface absorptions caused by water and con- 
tamination. Thus, in this section, we use r a b s = 329 ns 
to simulate a case with faster thermal generation. As dis- 
cussed above, here we assume that the Qj nt of the cavity 
is limited by absorption and not by losses to the outside of 
the cavity, since it appears to be the worst (but realistic 
case) for toroid microcavities. We employ r a b s , which is 
derived from the highest experimental Q m t ("Tabs = 329 ns 
corresponds to Q m t — 4 x 10 8 ). If we can clarify the re- 
quirements for demonstrating a Kerr bistable memory 
under this condition, it is a significant step toward the 
experimental realization of a Kerr bistable memory in 
silica toroid microcavities. 

First, in a similar way to that shown in Fig. (Ha), we 
employ a rectangular pulse to obtain the refractive index 
change AnKerr and Anxo- The calculation results are 
shown in Fig. [3] for three different r coup 2 values (t coup i 
is adjusted to satisfy T~l pl = T~l p2 + t^). Figure [3] 
shows that Atito is larger than Anxcrr in all three cases 
when t is larger than ~ 2.3 fj,s. This number gives us 
the upper limit of the Kerr memory holding time with- 
out the memory suffering from the TO effect. It also tells 
us that this number is insensitive to r coup 2. This result 
is consistent with Eq. (|5|) obtained from a simple model, 
where the equation is dependent on r a bs but independent 



of 



coup ■ 



Figure also shows the effect of the different 



charging speed resulted by different r coup 2 values. The 
cavity charging time is much faster for r coup 2 = Tj nt /100, 
which allows the cavity to reach a plateau AnKerr domain 
much faster. This enables us to have a longer "Kerr mem- 
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Fig. 2. (Color online) (a) AriKerr and Anxo vs. time when a 6 /xW rectangular pulse is inputted. The wavelength 
detuning S between the input light and the initial cavity resonance is 27 fm. T COU p2 is set equal to T lnt . (b) P ou ti 
vs. P in of a triangular input pulse with different 5 values, where P ou ti = |s ou ti| 2 - (c) P ou t2 vs. P; n with different 5, 
where P ou t2 = |s ou t2| 2 - 5 = 13, 20, and 27 fm correspond to (v / 3/2)AAfwhm, (3v / 3/4)AAfwhm, and V^AAfwhm, 
respectively. AAfwhm is the full-width at half-maximum of the cavity resonant spectrum width. (AAfwhm — 15 fm). 
(d) Optical bistable memory operation. The solid, dotted and broken lines represent p n , P ou ti and P ou t2, respectively. 



ory usable" regime, which allows us to use the cavity for 
longer as an optical Kerr memory. Figure [3] shows that we 
can maximize the Kerr dominant "Kerr memory usable" 
regime by setting r coup 2 as small as possible. Again, this 
result is in consistent with that obtained from the simple 
analysis described in § [2] 

Next, we show on how we can obtain Kerr optical bista- 
bility in a large absorption condition. Figure Ufa) and (b) 
show the input and output power relationships (P ou ti and 
P ou t2) for different r coup 2 values when a triangular pulse 
is inputted. As shown in Fig. SJa) and (b), a hysteresis 
loop is observed when r coup 2 = Tint/100, but the loops are 
deformed when r coup 2 > Ti nt /10. This is because the light 
cannot charge and discharge the cavity quickly enough 
before the heat accumulates in the system when r coup 2 
is large. Again, we obtained a clear hysteresis loop only 
when we made the coupling strong (i.e. T COU p2 is small). 

Finally, the memory operation for various r coup 2 values 



is shown in Fig. [SJ Since the response speed of P ou ti 
and P ou t2 depend on the total photon lifetime r to t, we 
normalized the temporal axis t by Ttot • Figure [5] shows 
clearly that the reset pulse does not work when r coup 2 = 
Ti n t and Tnt/10, and the memory operation cannot be 
obtained under this condition. If the system is operating 
in the Kerr dominant regime, we should be able to reset 
the state by injecting a negative reset pulse. The pulse 
width needed for the reset pulse is > T to t, since we can 
discharge the cavity within this time. However, Fig. [5] 
shows that significant heat is accumulating in the system, 
which prevents the system from resetting because Anxo 
cannot be reset by such a short negative pulse due to its 
much longer relaxation time. 

On the other hand, when T coup 2 < Tnt/100, we can suc- 
cessfully set and reset the system, and use the device as a 
Kerr bistable memory. However, the TO effect cannot be 
eliminated completely even in this case, and thus a hold- 
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Fig. 3. (Color online) AriKerr (solid line) and 
AriTO (broken line) vs. time for different 

7"coup2 

values. The input pulse power is set to 500 x 
(Oto° up2=Tint/100 /Qtot) 2 in fJ,W. The detuning S is set 
at \/3AAfwhm, where AAfwhm is 15 fm, 85 fm, and 
782 pm when T coup2 is r m t, Tmt/10, and T int /100, respec- 
tively. 



Fig. 4. (Color online) Input-output response of a trian- 
gular pulse input, (a) P ou ti vs. p n for different T coup 2. 
(b) P ou t2 vs. P; n for different T coup2 . The s-axis and y- 
axis are normalized as P™?f = P n u ?f = AT* = l - Thc 
detuning 5 is set at AAfwhm- The input pulse width 
is 15 fj,s, 3 fis and 0.3 tts and the pulse height is 4.5 iiW, 
160 /iW, 15 mW when T coup2 is r int , r int /10, and r int /100, 
respectively. 



ing time exists. P ou ti is automatically switched from ON 
to OFF, or vice versa for P ou t2, due to the thermal accu- 
mulation. Figure [S] shows that the memory holding time 
for the realistic case is about 500 ns. This value is in- 
consistent to the length of "Kerr memory usable" regime 
shown in Fig. [3] since the thermal nonlinearity depends 
on the accumulation of £/ a b s while the Kerr nonlinearity 
depend on the instantaneous U p . 

In this section we achieved Kerr bistable memory op- 
eration and obtained a sufficiently long memory holding 
time of 500 ns by allowing the system to charge and dis- 
charge quickly by adjusting r coup2 . 

5. Discussion: Power consumption 

We showed in previous sections that Kerr bistable oper- 
ation is possible by adopting a large coupling constant. 
The basic idea is to allow the light to charge and dis- 
charge before it turns into heat. However, stronger cou- 
pling with waveguides (i.e. smaller T COU p2) results in a 
lower total Q, which decreases the photon density in a 
cavity and makes the nonlinearity small. Hence, there is 



a trade-off between operating speed and operating power. 

Now we consider two measures for evaluating the loss 
of our system; Q lnt and U cons . Q m t is the figure of merit 
of a cavity, whose value gives the cavity loss, and U cons is 
the energy consumed for the operation. We often want to 
increase Qint and reduce U cons to build a lossless system. 
Table [T] compares these values for different systems. A 
bistable memory based on an ultrahigh-Q silica toroid 
microcavity with Kerr nonlinearity has a clear advantage 
over other schemes in terms of both losses. This lossless 
nature of this memory is the advantage of this system. 

First, a system with a high Q; nt yields a low loss, be- 
cause Qi n t corresponds to the fundamental loss charac- 
teristics of a cavity. As shown in Table [TJ the Qi n t of a 
toroid cavity is much higher than that of other types 
of cavities. Furthermore, we would like to note that 
both linear loss and nonlinear losses are significant in 
other systems, especially those systems that use TO and 
carrier-based nonlinearities to achieve bistability. This is 
because carrier generation is unavoidable. Free carrier 
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Table 1. Comparison of the power consumption of different systems 



System typo 


Material 


Nonlincarity 


Qint 


Pdrivc 








Exp./Cal. 


Rcfs. 


Microring cavity 


Si 


Thcrmo-optic 


1.43 x 10 5 


800 /iW 


3.7 x 


10" 


12 .7 


Exp. 


Ref. [11] 


Photonic crystal 


InGaAsP 


Carrier-plasma 


1.3 x 10 5 


250 /iW 








Exp. 


Ref. [13] 


Photonic crystal 


Semiconductor 


Kerr 


557 


133 mW 


2.0 x 


10" 


17(Mj 


Calc. 


Ref. [7J 


Toroid microcavity Si02 


Kerr 


4 x 10 8 


7.3 mW 


2.6 x 


10" 


23(6) j 


Calc. 


This work. 



a Estimated from the cavity resonance shift of (1/2)AAfwhm- We assume no TO and no carrier diffusion. 
b The energy change caused by the resonant wavelength shift AU = E/hwhm((1/2)AAfwhm/(-^0 + (1/2)AApwhm)) i s regarded as 
t^cons, where Chwhm is the energy that can cause a resonant wavelength shift of (l/2)AApwHM by using the Kerr effect. 
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Fig. 5. (Color online) Optical memory operation, when 
8 is \/3Afwhm • The input power of the drive light is 2 
/iW, 80 iiW, and 7.3 mW, when r coup2 is r int , r int /10 and 
Tint/ 100, respectively. The horizontal axis is normalized 
with the photon lifetime r to t- 



absorption (FCA) significantly decreases both the Qi nt 
and the transmittance of the system [36]. On the other 
hand, a Kerr bistable memory does not generate carriers 
and hence the system does not suffer from FCA loss. 

Secondly, U cons is extremely low in our system because 
Kerr nonlinearity does not absorb photons. U cons is the 
energy consumed for a bistable operation. We can es- 
timate [/ C ons from the required refractive index change, 
which is needed to obtain the cavity resonance shift for 
the operation. The values are large in memories based 
on TO and carrier-plasma effects even though they have 



large nonlinear efficiencies, which is because they incor- 
porate photon absorption. On the other hand, the total 
photon number for the Kerr effect remains unchanged af- 
ter the operation. Only the wavelength (energy) of the 
photons changes, which results in an extremely low en- 
ergy consumption. Hence U cons is almost negligible as 
shown in Table [TJ 

As shown in Fig. 5, we need P^" vc — 7.3 mW to drive 
the system, which is significantly smaller than the value 
shown in Ref. [7] , which uses the same Kerr nonlinearity. 
This is due to the high Qi nt of our system. However, this 
value is still larger than the experimentally demonstrated 
value in a photonic crystal nanocavity using carrier-based 
nonlinearity [hi} . Since the coefficients of carrier nonlin- 
carities in semiconductor materials are normally larger 
than the Kerr nonlinearity in Si02, it is not an easy 
task to reduce the driving power of our system to the 
same level. However, glass performs low propagation loss 
[37L Hl| , which is usually difficult to obtain with devices 
made of semiconductors. 

Although the driving power Pdrive of our system is not 
necessarily the smallest, the linear, nonlinear and con- 
sumption losses are significantly smaller than those of 
other devices. With this in mind, this device is attractive 
for applications that require high efficiency 77 = Pout /-Pin, 
such as quantum information processing [17]. 

6. Conclusion 

We rigorously modeled the Kerr and the TO effects in a 
silica toroid microcavity by combining CMT and FEM. 
We gained a clear understanding of the impact of adjust- 
ing the coupling, and showed that Kerr optical bistable 
memory operation is possible by adjusting the coupling 
between the cavity and the waveguides. The memory 
holding time was about 500 ns. Although the driving 
power was 7.3 mW, the energy consumed by the the sys- 
tem was extremely low. This is because unlike other 
nonlinearitics such as carrier of TO effects Kerr nonlin- 
earity does not absorb photons. In addition due to the 
ultrahigh-Q of the system, the energy loss outside the 
system is also low. Our Kerr bistable memory in a sil- 
ica toroid microcavity exhibits extremely low loss and 
thus is suitable for applications such as quantum signal 
processing. 
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